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ABSTRACT 

We study the relation between the rms mass fluctuations on 8/i _1 Mpc scales and f2 m using the recent 
clustering results of XMM-JVcwton soft (0.5-2 keV) X-ray sources, which have a median redshift of z ~ 1.2. 
The relation can be represented in the form erg = 0.34(±0.01)f2~ 7 where 7 = 7(O m , w) and it is valid for 
all w < —1/3 models. By combining the X-ray clustering and SNIa data we find that the model which 
best reproduces the observational data is that with: f2 m ~ 0.26, w ~ —0.90 and as ~ 0.73, which is in 
excellent agreement with the recent 3-year Wilkinson Microwave Anisotropy Probe results. 

Subject headings: cosmology:cosmological parameters - large-scale structure of universe 



1. INTRODUCTION 

The combination of the recently acquired, high quality, 
observational data on galaxy clustering, the SNIa Hub- 
ble relation and the CMB fluctuations, strongly support 
a universe with flat geometry and a currently accelerated 
expansion due to the combination of a low matter density 
and a dark energy component (eg. Riess, et al. 1998; Perl- 
mutter ct al. 1999; Percival et al. 2002; Efstathiou et al. 
2002; Spergel et al. 2003; Tonry et al. 2003; Schuecker et 
al. 2003; Riess et al. 2004; Tegmark et al. 2004; Seljak 
et al. 2004; Allen et al. 2004; Basilakos & Plionis 2005; 
Blake et al. 2006; Spergel et al. 2006; Wilson, Chen & 
Ratra 2006, for a review see also Lahav & Liddle 2006). 

From the theoretical point of view various candidates 
of the exotic "dark energy" have been proposed, most 
of them described by an equation of state pq = wpQ 
with w < —1/3 (see Peebles & Ratra 2003 and references 
therein). Note that a redshift dependence of w is also pos- 
sible but present measurements are not precise enough to 
allow meaningful constraints (eg. Dicus & Repko 2004; 
Wang & Mukherjee 2006). From the observational point 
of view and for a flat geometry, a variety of studies indi- 
cate that w < —0.8 (eg. Tonry et al. 2003; Riess et al. 
2004; Sanchez et al. 2006; Spergel et al. 2006; Wang & 
Mukherjee 2006 and references therein) 

Another important cosmological parameter is the nor- 
malization of the CDM power spectrum in the form of the 
rms density fluctuations in spheres of radius 8/i Mpe, 
the so called a%. A tight relation between a% and the 
f2 m has been derived mainly using the cluster abundance 
with a s ^ 0.52^ m 052 for a A cosmology (Eke, Cole & 
Frenk 1996). Also, Wang k, Steinhardt (1998) generaliz- 
ing to take into account dark energy models (with w > — 1) 
found: as =s 0.5Q-°- 21 +°- 22 ™-°- 33 ^. 

In this letter we use the clustering of high-z X-ray AGNs 
to estimate a new normalization of the CDM spectrum, 
valid for spatially flat cosmological models and also for 
w < — 1 (the so called Phantom models). Finally, com- 
bining our results with SNIa data (Tonry et al. 2003), we 
put strong constraints on the value of the equation of state 
parameter. 



2. X-RAY AGN CLUSTERING 

In a previous paper (Basilakos et al 2005) we derived the 
angular correlation function of the soft (0.5-2 keV) XMM 
X-ray sources using a shallow (2-10 ksec) wide-field sur- 
vey (~ 2.3 deg 2 ). A full description of the data reduc- 
tion, source detection and flux estimation are presented in 
Georgakakis et al. (2004). Here we describe only the basic 
points. The survey contains 432 point sources within an 
effective area of ~ 2.1 deg 2 (for f x > 2.7 x 10~ 14 erg cm~ 2 
s _1 ), while for f x > 8.8 x 10" 15 erg cm -2 s -1 the effective 
area of the survey is ~ 1.8 deg 2 . The details of the corre- 
lation function estimation, the various biases that should 
be taken into account (the amplification bias and integral 
constraint), the survey luminosity and selection functions 
as well as issues related to possible stellar contamination 
are presented in Basilakos et al. (2005). The redshift se- 
lection function of our X-ray sources was derived using 
the soft-band luminosity function of Miyaji, Hasinger & 
Schmidt (2000), and assuming the realistic luminosity de- 
pendent density evolution of X-ray AGNs and it predicts a 
characteristic depth of z ~ 1.2 for our sample (for details 
see Basilakos et al. 2005). 

Our aim here is to investigate the relation between the 
normalization of the CDM power spectrum (erg) and £l m 
in flat cosmologies with w < —1/3. A through study of the 
theoretical clustering predictions from different flat cosmo- 
logical models to the actual observed angular clustering of 
distant X-ray AGNs was presented in Basilakos & Plio- 
nis (2005). For the purpose of this study we use Limber's 
formula which relates the angular, w(0), and the spatial, 
£(r), correlation functions. In the case of a spatially flat 
Universe, Limber's equation can be written as: 



w(0) 



.Ho 
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E(z)dz / £(r,z)du (1) 



with e(z) = [n m (i + z) 3 + (i - n m )(i + 2 )3(i+-)]i/2. 

Also r is the physical separation between two sources, 
having an angular separation, 9, given by r ~ (1 + 



i 1/2 

z)^ 1 (u 2 + x 2 9 2 ) (small angle approximation). 
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number of objects within a shell (z, z + dz) and in a given 
survey of solid angle Q s is: 

^ = il s a?{z)n a <f>(x) (jL^j E -\z) . (2) 

where n s is the comoving number density at zero redshift 
and x(z) is the coordinate distance 

dy 



E(y) 



(3) 



Finally, the selection function <j){x) (the probability that 
a source at a distance x is detected in the survey) is esti- 
mated by integrating the appropriate Miyaji et al. (2000) 
luminosity function. 

2.1. The Evolution of Clustering 

It is well known (Kaiser 1984; Benson et al. 2000) that 
according to linear biasing the correlation function of the 
mass-tracer (£ t>j) and dark-matter one ((dm), are related 
by: 

£obj(r,z) = 6 2 (^DM(r,z) , (4) 

where b{z) is the bias evolution function. Here we use the 
bias evolution model of Basilakos & Plionis (2001; 2003), 
where we also compared in detail our model with that of 
Mo & White (1996) and Matarrese et al. (1997). As an ex- 
ample, our model predicts an AGN bias which is by <~ 30% 
greater and by ~ 35% lower than that of Mataresse et al. 
(1997) at z — and z = 3, respectively (see Basilakos et 
al. 2005). We remind the reader that our bias model is 
based on linear perturbation theory and the Friedmann- 
Lemaitre solutions of the cosmological field equations. For 
the case of a spatially flat cosmological model our general 
bias evolution can be written as: 

f°° (1 + v) 3 

b(z)=AE(z)+CE(z)J^ i^-dy + 1 . (5) 

Note that our model gives a family of bias curves, due 
to the fact that it has two unknowns (the integration 
constants A,C). The value of C is approximately found 
to be ~ 0.004, as was determined and tested in Basi- 
lakos & Plionis (2003). Note that E(0) = 1 and A = 

bo — 1 — C J °° gtfa) dt/, where b is the bias at the present 
time. We have tested the robustness of our results by in- 
creasing C by a factor of 10 and 100 to find differences of 
only ~ 5% in the fitted values of Q m and b . This behavior 
can be explained from the fact that the dominant term in 
the right hand side of eq. (5) is the first term [oc (l + z) 3 / 2 ] 
while the second term has a slower dependence on redshift 

[oc(l + *)]■. 

We quantify the underlying matter distribution cluster- 
ing by presenting the spatial correlation function of the 
mass £r>M{r,z) as the Fourier transform of the spatial 
power spectrum P(k): 



£r>M(r,z) = 



(1 + *) 



"(3+0 



where k is the comoving wavenumber. Note that the pa- 
rameter e parametrizes the type of clustering evolution 
(eg. de Zotti et al. 1990). In this work we utilize a clus- 
tering behavior which is constant in comoving coordinates 
(6 = -1.2). 

4 Likelihoods arc normalized to their maximum values. 



As for the power spectrum, we consider that of CDM 
models, where P(k) — Pok n T 2 (k) with scale-invariant 
(n = 1) primeval inflationary fluctuations (we verified that 
a small change of n, eg. n ~ 0.95 according to the 3-year 
WMAP does not produce appreciable differences on our 
results). In particular, we use the transfer function pa- 
rameterization as in Bardeen et al. (1986), with the cor- 
rections given approximately by Sugiyama (1995) while 
the normalization of the power spectrum is given by: 



Po 



2tt 2 oJ 



f 

Jo 



T\k)k n+z W\kR)dk 



(7) 



where as is the rms mass fluctuation on R — 8h x Mpc 
scales and W(kR) is the window function given by 



W(kR) = 



3(sinfci? — kRcoskR) 



(8) 



(fci?) 3 

Note that we also use the non-linear corrections introduced 
by Peacock & Dodds (1994). 

3. COSMOLOGICAL CONSTRAINTS 

3.1. X-ray AGN Clustering likelihood 

Following the same notations as in Basilakos & Plionis 
(2005) in order to constrain the cosmological parameters 
we use a standard x 2 likelihood procedure to compare the 
measured XMM soft source angular correlation function 
(Basilakos et al. 2005) with the prediction of different spa- 
tially flat cosmological models. The likelihood estimator 4 
is defined as: 



£ AGN (c)cxexp[- X i GN (c)/2] 



with: 



Xagn 



w t h(0i,c) - W obs (0i) 



(9) 



where c is a vector containing the cosmological parame- 
ters that we want to fit and <7j is the uncertainty of the 
observed angular correlation function. We make clear that 
we work within the framework of a flat cosmology with 
primordial adiabatic fluctuations and baryonic density of 
Q h h 2 = 0.022(±0.002) (eg. Kirkman et al. 2003; Spergel 
et al. 2006). Also utilizing the results of the HST key 
project (Freedman et al. 2001) we fix the Hubble constant 
to its nominal value; ie., h ~ 0.72, derived also by our pre- 
vious AGN clustering analysis (Basilakos & Plionis 2005). 
In that work the la uncertainty of the marginalized value 
of h was found to be only <~ 0.03. Note that since we fix 
in the following analysis the values of both h and fib, we 
do not take into account their quite small uncertainties. 

The corresponding vector that we have to fit is c = 
(fi m , w, as, bo). In this paper, we sample the various pa- 
rameters as follows: the matter density f2 m G [0.01,1] 
in steps of 0.01; the equation of state parameter w e 
[—2, —0.35] in steps of 0.05, the rms matter fluctuations 
as G [0.4, 1.4] in steps of 0.02 and the X-ray sources bias 
at the present time 6 G [0.5, 3] in steps of 0.05. Note that 
in order to investigate possible equations of state, we have 
allowed the parameter w to take values below -1. Such 
models correspond to the so called phantom cosmologies 
(eg. Caldwell 2002; Corasaniti et al. 2004). 



Constraining the CDM spectrum normalization 
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The resulting best fit parameters are presented in Ta- 
ble 1. It is important to note that our estimate for the 
as parameter is in very good agreement with that derived 
(er 8 ~ 0.75) by the recent 3-years WMAP results (Spergel 
et al. 2006). In Fig.l we present the la, 2a and 3a confi- 
dence levels in the (Q m ,a$), (as, bo) and (fl m ,bo), planes 
by marginalizing the first one over b and w, the second 
one over £l m and w and the latter over as and w; while 
the dot in Fig. 1 corresponds to the best fitted values 5 . 




0.2 0.4 0.6 0.B 



Fig. 1. — Likelihood contours in the following planes: (f2 m ,crg) 
(upper right panel), (<rg,&o) (upper left panel), (f2 m ,6o) (bottom 
left panel) and the (S7 m ,w) (bottom right panel). The contours are 
plotted where — 21n£/£ max is equal to 2.30, 6.16 and 11.83, cor- 
responding to lcr, 2a and 3<r confidence level. Finally, the thick 
contours corresponds to the SNIa likelihoods. 

Therefore, utilizing the clustering properties of our 
XMM sources and allowing for the first time values w < — 1 
(Phantom models) we derive the (^l m ,as) relation, which 
can be fit (within the la uncertainty) by 



as = 0.34(±0.01) fl m ^ n " 



,w) 



(10) 



with 

7(O m ,w) = 0.22(±0.04)-0.40(±0.05)io-0.052(±0.040)fi m . 

In Figure 2 we present the results of the likelihood anal- 
ysis for different values of w (points with errorbars) and 
the previous fit as a continuous line. 

Note that eq. (10) produces as values which are signifi- 
cantly smaller than the usual cluster normalization (Wang 
& Steinhardt 1998) but are in good agreement with the 
3-years WMAP results; for example for w ~ —1 and 
O m ~ 0.28 we get a s =i 0.73±0.03. It should be mentioned 
that in our previous work (Basilakos & Plionis 2005) we 
had imposed a high as normalization, based on the cluster 
abundance, while here we leave the as parameter free. 

The lower right panel of Fig. 1 shows the la, 2a and 3cr 
confidence levels (continuous lines) in the (fi m , w) plane by 
marginalizing over the as and the bias factor at the present 
time. It is evident that w is degenerate with respect to fi m 

5 Hereafter, when we marginalize over the equation of state parameter 



and that all the values in the interval —2 < w < —0.35 are 
acceptable within the la uncertainty. Therefore, in or- 
der to put further constraints on w we additionally use a 
sample of 172 supernovae (see Tonry et al. 2003). 

3.2. The AGN+SNIa likelihoods 

Here we combine the X-ray AGN clustering properties 
with the SNIa data by performing a joined likelihood anal- 
ysis and marginalizing the X-ray clustering results over as 
and bo (see Table 1) and thus the vector c now becomes: 
c = (fl m ,w). The SNIa likelihood function can be written 



as: 



£ bNla (c)«e X ph X I NIa (c)/2] 



with: 

xlNia(c) = 



172 r 



logl^te.cJ-logD^te) 



(11) 



where Dj^(z) is the dimcnsionlcss luminosity distance, 
D-l(z) = H Q (1 + z)x(z) and z% is the observed redshift. 
The thick lines in Fig. 1 represents the la, 2a, and 3a, 
confidence levels in the (Q m ,w) plane. We find that the 
best fit solution is f2 m = 0.30 ± 0.04 for w > — 1, in com- 
plete agreement with previous SNIa studies (Tonry et al. 
2003; Riess et al. 2004). 

We now join the likelihoods: 

£ joint (f! m , w ) = £ AGN x£ SNIa , 

which peaks at fl m = 0.26±0.04 with w = -0.90^;^. Us- 
ing eq. (10) we find that the normalization of the power 
spectrum that corresponds to these cosmological param- 
eters is as — 0.73. It should be pointed out that our 
results are in excellent agreement with those derived by 
Spergel et al. (2006) using the recent WMAP (3-years) 
data: n m ~ 0.24, w ~ -0.97 and a s ~ 0.74. 




Fig. 2. — The (n m ,(78) plane using different values for the equa- 
tion of state parameter (points). The errors correspond to lcr (2.30) 
confidence levels. The continuous line corresponds to the fit given 
by equation 10. 

we will use w = — 1. 
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Many other recent analyzes utilizing different combina- 
tions of data seem to agree with the former cosmologi- 
cal model. For example, Sanchez et al. (2006) used the 
WMAP (1-year) CMB anisotropics in combination with 
the 2dFGRS power spectrum and found f2 m ~ 0.24 and 
w ~ —0.85, while Wang & Mukherjee (2006) utilizing the 
3-years WMAP data together with SNIa and galaxy clus- 
tering results found w ~ —0.9 (see also Wilson et al. 2006). 

4. CONCLUSIONS 

We have combined the clustering properties of distant 
X-ray AGNs, identified as soft (0.5-2 keV) point sources 
in a shallow ~ 2.3 deg 2 XMM survey, with the SNIa data. 



From the X-ray AGN clustering likelihood analysis alone 
we have estimated the normalization of the CDM power 
spectrum and find that the rms density fluctuation in 
spheres of radius 8h~ 1 Mpc is fitted by: 

<j s ~ 0.34(±0.01)f7 m - 22 + - 40 ™+°- 052O » 

which is valid also for Phantom models (w < —1). Fur- 
thermore, a joined likelihood analysis between the X-ray 
and SNIa data provides a best model fit with: £l m ~ 0.26 
and w ~ —0.90, which corresponds to er 8 ~ 0.73, in agree- 
ment with the recent 3-years WMAP results (Spergel et 
al. 2006). 
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Table 1 

cosmological parameters from the likelihood analysis: the 1 st column indicates the data used (the last row 
corresponds to the joint likelihood analysis). errors of the fitted parameters represent la uncertainties. 
Note THAT FOR THE JOINED ANALYSIS THE CORRESPONDING results are marginalized OVER THE <7 8 and the bias FACTOR 

AT THE PRESENT TIME, FOR WHICH WE USE THE VALUES INDICATED. 
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